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Abstract. Coccolithophore thanatocoenosis from 33 surface sediment samples in the Argentinean Continental Margin were analyzed. Emiliania 
huxleyi, the most ubiquitous species in today’s ocean, occurs at relative abundances higher than 60% and tends to obscure the variations in 
abundance of less-represented species. This study aims to compare two intensively used counting methods in coccolithophore studies, ie, 
including and excluding E. huxleyi from the data set, and thus compare which of both methodologies better reflects the known oceanographic 
conditions in the area. In the two data sets, species abundances were converted to percentages and Cluster Analysis were performed. The 
result including £. huxleyi reveals two groups (A and B) separated by compositional Euclidean distance of 12.8. Group A is restricted to the 
southeastern area (near Malvinas Islands) and is dominated by Gephyrocapsa muellerae (> 56%). Group B is dominated by E. huxleyi (> 93%) and 
is distributed into two areas: in the southwest (near Tierra del Fuego) and the north (between 40 °S—47 °S). The result excluding E. huxleyialso 
shows two groups (A’ and B’) distant to a compositional Euclidean distance of 15.8. Group A’ has both a southern and northern distribution, 
whereas Group B’ is present only in the northern area. These two groups cannot be explained by the relation between the coccolithophore taxa 
dominance and the surface oceanographic conditions. Our results indicate that the inclusion of E. hux/eyiis the most appropriate approach 
since it seems to better correlate with the known oceanographic regime in the Argentinean Continental Margin. 


Key words. Calcareous nannoplankton. Counting methodology. Data analysis. Surface sediment. Southwest Atlantic Ocean. 


Resumen. EVALUACION DEL MEJOR CONJUNTO DE DATOS PARA EL ANALISIS ESTADISTICO DE LA DISTRIBUCION DE COCOLITOFORIDOS EN 
EL MARGEN CONTINENTAL ARGENTINO. Se analizé la tanatocenosis de cocolitoféridos en 33 muestras de sedimentos superficiales del Mar- 
gen Continental Argentino. Emiliania huxleyi, especie actual mas ubicua, se encuentra en abundancias relativas superiores al 60% y tiende a 
ocultar las variaciones de especies menos representadas. El! objetivo de este estudio es comparar dos métodos de conteo de uso intensivo, ie., 
incluir y excluir E. huxleyi del conjunto de datos y asi comparar cual de ambas metodologias refleja mejor las condiciones oceanograficas. En los 
dos conjuntos de datos, las abundancias de los taxones se convirtieron en porcentajes y se realizaron Analisis de Agrupamiento. El resultado 
que incluye E£. huxleyi revela dos grupos (A y B) separados por una distancia euclidiana de 12,8. El grupo A, restringido a la zona sureste (cerca 
de las islas Malvinas), esta dominado por Gephyrocapsa muellerae (> 56%). El grupo B, dominado por E. huxleyi(> 93%), se distribuye en dos areas: 
suroeste (cerca de Tierra del Fuego) y norte (entre 40 °S—47 °S). El resultado que excluye a £. huxleyi muestra dos grupos (A’ y B’) auna dis- 
tancia euclidiana de 15,8. El grupo A’ tiene una distribuci6n sur y norte, mientras que el grupo B’ esta presente en el norte. Estos dos agrupa- 
mientos no pueden explicarse por la relacion entre la dominancia de los cocolitoféridos y las condiciones oceanograficas superficiales. Nuestros 
resultados indican que la inclusi6n de £. huxleyi es el enfoque mas apropiado ya que parece correlacionarse mejor con el régimen oceanogra- 
fico conocido en el Margen Continental Argentino. 


Palabras clave. Nanoplancton calcareo. Metodologia de conteo. Analisis de datos. Sedimentos superficiales. Océano Atlantico Sudoeste. 


CoccoLitHoPHoRES, phytoplankton belonging to division — of calcified scales, the coccoliths, which form a major com- 
Haptophyta and class Prymnesiophyceae, are one of the ponent of pelagic carbonate oozes and are abundant in 
most important primary producers of the oceans (Young eta/, | sea-floor sediments above calcite compensation depth. 


2003). A distinctive feature of these algae is the presence _ Isolated coccoliths can preserve the major compositional 
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characteristics of the overlying photic-zone communities 
due to protected and accelerated sedimentation inside fecal 
pellets or marine snow (Steinmetz, 1994). Coccolithophores 
are present in a wide range of marine environments ranging 
from tropical to subpolar latitudes. Biogeography of living 
coccolithophores responds sensitively to changes in 
environmental parameters but there is still much discussion 
as to which conditions control the vertical and horizontal 
distributions (Winter et a/, 1994). However, their diversity 
and distribution in the photic zone seem to be affected by 
surface oceanic circulation and parameters such as sea 
surface temperature, salinity and nutrient availability 
(Winter & Siesser, 1994; Ziveri et a/, 2004). Based on those 
observations, coccolith assemblages in surface sediment 
have been used as proxies to monitor variations in 
paleoceanographic conditions (e.g., Brand, 1994; Winter & 
Siesser, 1994; Saavedra-Pellitero et a/, 2013). 

Knowledge of the spatial distribution of living coccol- 
ithophores is crucial to understand their biology and envi- 
ronmental preferences. Significant biogeographic contri- 
butions of modern coccolithophores in the Atlantic Ocean 
have been made by Mcintyre & Bé (1967) and Okada & 
McIntyre (1977, 1979). McIntyre & Bé (1967) established a 
latitudinal biogeographical nannoplankton zonation based 
on temperature ranges and assemblage differences. The 
four floral zones (tropical, subtropical, transitional, and sub- 
antarctic or subarctic) are typified in the warmest and 
coldest by the dominance of one or more key species. The 
transitional flora, which is localized in the ACM, is charac- 
terized by the dominance of placoliths such as Gephyrocapsa 
ericsonii, Emiliania huxleyi and Calcidiscus leptoporus, and 
lesser amounts of Coccolithus pelagicus. However, these 
latitudinal zonations do not consider the influence of coastal 
currents, eddies and upwelling, concealing species-specific 
biogeography. Mostajo (1986) analyzed coccoliths in sedi- 
ment samples from the Southwest Atlantic Ocean including 
the ACM and identified the Subtropical, Transitional and 
Subantarctic water masses sensu Boltovskoy (1970). Be- 
sides, she defined the boundaries of the Confluence Zone (= 
Convergence) for summer and winter. More recently, 
Boeckel et al. (2006) analyzed the coccolith composition of 
sediment samples from the South Atlantic and the Southern 


Ocean concerning to different environmental parameters. 
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According to them, coccolith distribution is related to the 
nutricline and thermocline position. 

In most studies dealing with the distribution of modern 
coccolithophores, Emiliania huxleyi was found to be the 
dominant species (e.g, Okada & McIntyre, 1979; Baumann et 
al, 2008; Saavedra-Pellitero & Baumann, 2015). In addition, 
Emiliania huxleyi has a wide biogeographic range, being the 
most ubiquitous coccolithophore in today’s oceans. In 
subarctic and subantarctic waters £. huxleyi constitutes 
almost 100% of the coccolithophore assemblage (Mcintyre 
& Bé, 1967) and tends to overlap the relative abundance 
fluctuations of other less common taxa. As a result, for 
ecological interpretations and paleoceanographic recons- 
tructions, some studies chose to remove E. huxleyi from 
quantitative analysis (e.g, Giraudeau, 1992; Giraudeau & 
Rogers, 1994; Boeckel et a/., 2006), whereas other studies 
decided to include it (eg, Boeckel & Baumann, 2004; 
Saavedra-Pellitero et a/, 2013). However, whereas the 
exclusion of £. huxleyi seems to be a smart decision to 
highlight distribution patterns from other less represented 
but stenotopic species in the assemblages, the impact of its 
exclusion from the data sets had never been tested. This 
situation encouraged us to test which data set better 
correlates with known oceanographic conditions in the 
ACM, prior to focusing our studies on the environmental 
parameters controlling its distribution and reconstructing 
paleoceanographic conditions in the fossil record. 

The main goal of this study is to compare two different 
counting methodologies —including and excluding E. hux/eyi- 
by defining two data sets from a same sample set from the 
ACM (40 °S—55 °S and 66 °W-55 °W) and explain which of 
both methodologies better reflects the known oceano- 
graphic conditions. This study is an important prerequisite 
for future paleoecologic and paleoceanographic studies in 
sediments from the ACM, since it will allow the comparison 
between fossil (Pleistocene—Holocene) and recent assem- 


blages with a high degree of confidence. 


OCEANOGRAPHIC SETTING 

The South Atlantic Ocean plays a particular role in the 
thermohaline circulation and water mass distribution, 
making it an important region of nutrient exchange. The 


South Atlantic Ocean is dominated by three water zones 
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Figure 1. Upper ocean circulation and Sea Surface Temperature (SST in °C, Boyer et a/, 2013) annual average at O m in the southwestern 
Atlantic Ocean. Currents are indicated with a white line, according to Piola & Matano (2017) and Hernandez-Molina et a/. (2010), and isotherms 
with a black line. Oceanographic fronts are indicated with a dash line according to Peterson & Stramma (1991). Legend for the physiographic 
reference points, in alphabetical order: BB - Burdwood Bank, MI - the Malvinas Islands, TF - Tierra del Fuego. Plotted with Ocean Data View 
(ODV) software version 5.1.7 (Schlitzer, 2015). 
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divided by three oceanographic fronts: Subtropical Front, 
Subantarctic Front and Antarctic Polar Front (Peterson & 
Stramma, 1991) (Fig. 1). The Subtropical Zone is distin- 
guished by warm, salty and poor-nutrient waters. To the 
south is the Subantarctic Zone, where the ACM is located, 
which is characterized by an abrupt decline in salinity and 
temperature of surface waters. The Antarctic Polar Zone is 
dominated by waters with very high nutrient content and 
sea surface temperature lower than 10 °C. 

In the Southwestern Atlantic Ocean, the ACM is a highly 
significant region in the oceanographic-climatic system 
(Violante et a/, 2017). The ACM, one of the largest conti- 
nental margins on Earth, has different environments (shelf, 
slope, submarine canyons and rise) and the oceanic circula- 
tion is characterized by the opposing flows and confluence 
between the Malvinas and the Brazil currents (Matano et al, 
2010) (Fig. 1). The MC is a narrow branch of the ACC witha 
northward flow along the continental slope of Argentina 
near 34-38 °S. It is characterized by cold (T< 15 °C) and 
“fresh” waters (S< 34.2 PSU). This current is highly ho- 
mogenous and barotropic, thus it has less thermohaline 
stratification than waters of tropical or subtropical origin 
(Matano, 1993; Piola & Matano, 2017; Violante et a/, 2017). 
Most of the ACC water flowing eastward in the northern 
Drake Passage loops northward to form the MC. A portion of 
the upper flow deflects northward west of the Burdwood 
Bank and the remainder of the northern ACC deflects north- 
ward across a gap located east of the Burdwood Bank. Both 
branches of the MC rejoin to the north of the bank (Piola & 
Matano, 2017). In the northern portion of the Drake Passage, 
the mean annual surface temperature of the MC is around 4 
°C and increases northward up to 16°C where the currents 
separate from the continental boundary. The MC percolates 
over the shelf influencing the regional ecosystems (Piola et 
al, 2010). At 38°S, this current collides with the tropical- 
warm (T> 20 °C) and salty waters (S> 36 PSU) of the BC 
flowing towards the south along the Brazilian continental 
margin, originating a highly dynamic frontal zone: the BMC 
(Gordon, 1989; Piola & Matano, 2017). The encounter of MC 
and BC generates a sharp horizontal and vertical gradient in 
temperature, salinity, density and nutrient content (Gordon, 
1989; Peterson & Stramma, 1991; Bianchi et a/., 1993; Piola 
& Matano, 2001; Matano et a/,, 2010). As a result of in- 
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creased upper-water column stability and nutrient avail- 
ability, the BMC is a region of high primary productivity and 
major sink for atmospheric CO, (Feely et a/,, 2001). The BMC 
shows spatial and seasonality fluctuation characterized by 
an equatorward (poleward) displacement of the front during 
austral winter (summer) (Boltovskoy, 1981; Mostajo, 1986). 


MATERIALS AND METHODS 

During three expeditions on board the A.R.A. Austral 
oceanographic vessel, from September 2017 to October 
2018; 33 sea-floor sediment samples were obtained from 
the ACM. According to the locations, stations were grouped 
in three geographic areas: north area (40 °S to 47 °S and 55 
°W to 60 °W), southeast area (53 °S to 54 °S and 61 °Wto 
64 °W, near the Malvinas Islands), and southwest area (55 
°S and 65 °W, near Tierra del Fuego) (Fig. 2, Tab. 2). The 
samples were obtained from gravity core or box corer and 
the uppermost centimeter of comparatively undisturbed 
sediment column was sampled. 

AMS “C dates were performed on six samples from 
gravity and box cores, with the intention of having a better 
age control of the three different station areas investigated. 
Dating was conducted on picked planktonic foraminifera 
Globigerina bulloides (d'Orbigny, 1826). The analysis was 
carried out at DirectAMS Radiocarbon Dating Service, 
Washington, USA. AMS “C ages were converted to calendar 
years using CALIB 7.8.2 software (Stuiver & Reimer, 1993), 
Marine20 calibration curve (Heaton et a/, 2020) and no AR. 
In addition, age and sedimentation rates from other sites 
were estimated from the adjacent core locations. Numerous 
box and multicores in the vicinity of the samples positions 
exhibit relatively continuous late Quaternary sedimentation 
rates in which top sediments are allocated to Holocene 
(based on AMS “4C dates on planktic foraminifera, total 
organic carbon and alkenones, e.g., Benthien & Muller, 2000; 
Vink, 2004; Mollenhauer et a/., 2006; Ruhlemann & Butzin, 
2006). From the apparently coeval deposition of these 
sediments, we infer that transport processes occurred due 
to rapid advection in the upper water column (Benthien & 
Muller, 2000; Mollenhauer et a/, 2006). Sedimentation 
rates of ~13 cm per 1,000 years can be inferred for the 
Holocene hemipelagic sediments at a water depth of ~2,000 
m (Mollenhauer et a/., 2006). 
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Figure 2. Map of Argentinean Continental Margin with site number locations. The three station areas are indicated with rectangles: north area 
(dash line), southeast area (continuous line) and southwest area (dotted line). Exact positions are indicated in Table 2. Plotted with GMRT data 
(Ryan et al, 2009). 
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Sample preparation and coccolithophore counting 

The 33 samples were prepared for study under a 
polarized light microscope using a modification of the 
Standard Gravity Settling Technique of Bramlette & Sullivan 
(1961). In each case, 0.10 to 0.12 g of wet sediment was 
mixed in 10 ml of mineral water with the addition of 
hydrogen peroxide and hexametaphosphate. This solution 
was ultrasonically treated for 30 seconds to facilitate 
disaggregation of coccoliths from the sediment. After 30 
seconds, 200 ul of the suspension was mounted on a cover 
slide and affixed with Norland 61 optical adhesive. Coccolith 
observation was done using a Leica DM750 P microscope 
at 1,000x magnification. The preservational state of the 
coccoliths was examined under microscope. Samples are 
deposited and curated in the Y-TEC Laboratory of 
Biostratigraphy Micropaleontological Repository under the 
acronym YT.RMP_N. and numbers 000006.1 to 000006.35. 

For the two data sets (including and excluding E. huxleyi) 
at least 300 coccoliths were counted on an arbitrarily 
chosen transect, using the same sample set. This counting 
effort represents a 95% confidence of recording any species 
with a relative abundance of 1% (Thierstein et al, 1977), 
ensuring that most of the species were included (exceeds 
99% sensu Fatela & Taborda, 2002). Absolute abundance 
(coccoliths per field of view) was calculated for each sample. 
For a qualitative control, in the data set including E. huxleyi, 
after reaching 300 specimens, additional 150 fields of view 
were examined to identify any species present in the 
sample. Those species were not included in the Cluster 
Analysis. 

Coccoliths of Gephyrocapsa, which are smaller than 3 um 
and hard to identify due to their tiny size, were grouped 
under the name of Gephyrocapsa spp. small. Systematic 
criteria in this study follow Hine & Weaver (1998), Young et 
al, (2003), Jordan et al. (2004) and the online Nannotax 


Catalog (www.nannotax.org). 


Coccolith dissolution index and Shannon diversity 
index 

Carbonate calcium dissolution in deep waters can modify 
the taxonomic composition of recent coccolithophore assem- 
blages making difficult to conduct biogeographic and paleoe- 
cological reconstructions (Thierstein, 1980). To estimate the 
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effect of carbonate dissolution on the coccolith assemblages, 
the preservation index, CEX', from Boeckel & Baumann 
(2004) was applied. This index is based on the CEX from 
Dittert et al (1999), namely Calcidiscus leptoporus - Emiliania 
huxleyi Dissolution Index. The CEX’ index (Boeckel & Baumann, 
2004) is a ratio that compares the small and delicate pla- 
coliths of £. huxleyi and Gephyrocapsa ericsonii to strongly 
calcified coccoliths such as C. /eptoporus: 


CEX’=[E. huxleyi (%)+G. ericsonii (%)]/[ E. huxleyi (%) 


+G. ericsonii (‘)+C. leptoporus (%)] 


Carbonate dissolution has more effect on small 
placoliths than on robust ones and hence, the ratio of these 
taxa will change with increasing dissolution. According to 
Boeckel & Baumann (2004), a CEX’ value lower than 0.6 
indicates a depositional environment below the calcite 
lysocline. In this study, Gephyrocapsa spp. small (< 3 tm) 
was considered instead of G. ericsonii. 

H’ was used to measure the species diversity in each 
sample. This value varies from 0, when the communities 
have a single taxon, to higher values when communities 
have more taxa with an equitable number of individuals. H’ 
was performed with Past 4.02 (Hammer et a/., 2001). 


Statistical methods 

Coccolith assemblages are composed of species that 
appear frequently (common species), and others with very 
low frequency (rare species). Clifford & Stephenson (1975) 
considered that the probability of occurrence of a rare 
species in a sample is very small and therefore the co- 
occurrences lack of a spatial pattern. Thus, it is not possible 
to link them to habitat requirements because they occur 
randomly. Owing to these features, before the statistical 
analysis, species with an abundance less than 1% of the 
total specimens in all samples were considered rare species 
and were eliminated in both data sets, focusing on the 
variations of the most abundant taxa. 

For the statistical treatment of both data sets, 
coccolithophore relative abundances were expressed as 
percentages, as are known in Aitchison’s composition 
(Aitchison, 1986). The percentages of composition are 


relative values adding up to 100; thus, there is a statistical 
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dependence and spurious correlations between percen- 
tages. Aitchison (1986, 1999) proposed log-ratios strategy 
to break the dependences. Application of this strategy in 
coccolithophore studies can be found in Monechi et ai. 
(2000) and Buccianti & Esposito (2004). 

In order to compare the two counting methodologies, 
two hierarchical Cluster Analysis were performed. One 
widely used clustering algorithm is Ward’s method, which 
attempts to generate clusters to minimize the within- 
cluster variance. For each cluster, the means for all the 
variables are computed. Next, for each object, the squared 
Euclidean distance to the cluster means is calculated. 
These distances are summed for all the objects. This 
variance is a weighted squared distance between cluster 
centre. However, Euclidean distance is not appropriate for 
compositional data because it is not perturbation invariant 
and not sub-compositionally dominant. Therefore, to adapt 
Ward's method, data matrices were re-expressed using the 
clr as proposed by Martin-Fernandez et al. (1998) and 
Aitchison et al. (2000). 

The ci/ris isometric and was calculated with the following 


equation: 


clr(x) = (Inx1/g(x),In x2/g(X),... ; a 
In xd/g(x))=(z1,Z2,...,2d), where g(x) -[T 


where d, is the component number in the composition, 
clr (x)and zare the raw percentages and transformed value, 
respectively, for samples in rows and species in column 
matrix data, with the geometric mean g/x) of all components in 


the composition serving as the common divisor for each row. 


Null values were replaced using multiplicative replace- 
ment strategy with 0.1% (Martin-Fernandez et a/., 2003) 
before transformation to avoid division by 0, as recom- 
mended by Aitchison (1986). Clusters were identified with 
the geometric mean of species distribution after the closure 
operation was applied (the center in compositional words). 
Besides, the c/r variances of each species within the group 
were computed. 

Cluster Analysis and c/rtransformation were performed 
with Past 4.02 (Hammer et a/., 2001). 

Institutional abbreviations. A.R.A., Austral, Armada de la 
Republica Argentina Austral, Argentina; Y-TEC, Yacimien- 
tos Petroliferos Fiscales Tecnologia, Berisso, Argentina; 
YT.RMP_N., Yacimientos Petroliferos Fiscales Tecnologia, 
Repositorio Micropaleontolégico, Nanofdésiles. 

Oceanic abbreviations. ACM, Argentinean Continental 
Margin; ACC, Antarctic Circumpolar Current; BC, Brazil 
Current; BMC, Brazil-Malvinas Confluence; MC, Malvinas 
Current. 

Indices, statistical and dating abbreviations. AMS ™“C, Ac- 
celerator mass spectrometry radiocarbon; BP, before pres- 
ent; CEX and CEX’, Calcidiscus leptoporus - Emiliania huxleyi 
Dissolution Index; c/r, centered log-ratio transformation; H’, 
Shannon diversity index. 


RESULTS 
Radiocarbon dating 

The six calibrated AMS "“C ages obtained from sea-floor 
samples indicate Late Holocene/Recent sediments for the 
three investigated station areas (Tab. 1). Ages are given in 


™C radiocarbon and calendar years BP. 


TABLE 1 — Accelerator mass spectrometry radiocarbon dates and calibrated ages obtained on six sea-floor sediment samples. 


Site number Sample Direct AMS Code ™4C Radiocarbon age + 1 o error Calendar age BP + 1 o error 

3 Geo0_GC8 D-AMS 044169 811 ae 23} 124-428 276 + 152 
8 Geo01_BC11 D-AMS 044170 806 +21 120-425 272.5 + 152.5 
133 Geo01_GC16 D-AMS 044171 978 ee {I} 300-528 4142114 
21 Geo02_GC22 D-AMS 044172 2729 +23 2110-2412 2261+ 151 
14 Geo02_GC24 D-AMS 044173 2848 22 D2 2286-2614 2450 + 164 
25 Geo02_GC30 D-AMS 044174 1983 +21 1268-1512 1390 + 122 
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RIVAS Ef AL: COCCOLITHOPHORES COUNTING METHODOLOGY FROM ARGENTINA 


Coccolithophore assemblages in the ACM 

All the samples contain coccoliths and a total of 21 taxa 
were identified in both data sets (Tabs. 2, 3). Taxa are listed 
alphabetically with full citations in the Taxonomic Appendix. 
The most abundant coccolith taxa are F. huxleyi (77.54% 


average), Gephyrocapsa muellerae (16.53% average), Calcidiscus 


leptoporus (2.65% average) and Gephyrocapsa spp. small 
(1.58% average) (Fig. 3). Pre-Quaternary taxa, which 
occurred sporadically (Chiasmolithus spp., Helicosphaera sellii, 
Reticulofenestra pseudoumbilicus, Sphenolithus moriformis, 
Watznaueria spp.), were found in 22 samples and grouped 
as ‘Reworked’. Absolute abundance (coccoliths per field of 


TABLE 3 — Coccolith abundance chart, excluding Emiliania huxleyi. 
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Figure 3. Coccolithophore species considered for the Cluster Analysis with both data sets. 1, Emiliania huxleyi, sample YT.RMP_N.OO0006. 12; 
2, Gephyrocapsa muellerae, sample YT.RMP_N.000006.13; 3, Gephyrocapsa oceanica, sample YT.RMP_N.000006.16; 4, Gephyrocapsa spp. 
small, sample YT.RMP_N.000006.19; 5, Helicosphaera carteri, sample YT.RMP_N.000006.28; 6, Helicosphaera princei, sample 
YT.RMP_N.000006.11; 7, Umbellosphaera tenuis, sample YT.RMP_N.000006.23; 8, Calcidiscus leptoporus, sample YT.RMP_N.O00006.12. 


Scale bar: 10 pm. 


view) and richness showed the same trend: highest values 
in the northern sites and decreasing towards the south. H’ 
values range between 0.13 and 1.05 with highest values 
at the southern sites and lowest values at the northern 
sites. 

In all the samples, coccolith preservation is good. High 
CEX’ index values, all over 0.6 (min= 0.88, max= 0.99), 
indicate that dissolution is not significant. In almost all the 
samples, E. huxleyi is the dominant species except for a 
group located in the southeast area (which is dominated by 
G. muellerae). The ubiquitous E£. huxleyi exhibits the highest 
abundance in the northern sites (> 90%) and minimum in the 
southeastern sites (= 15%). 

Accordingly, maximum percentages of G. muellerae are 
recorded at the southeastern sites, close to the Malvinas 
Islands (60%-—70%). The distribution of Gephyrocapsa spp. 
small basically concurs with the G. muellerae trend. It shows 
the highest percentages (mean 3.4%) at the southern sites 
and, conversely, the minimum values (< 1%) at the northern 
sites. Calcidiscus leptoporus is especially abundant at the 
southwestern sites, near Tierra del Fuego (mean 7%), and 
less abundant at the southeastern and northern sites (= 2%). 
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Cluster Analysis 

Clusters Analysis were based on taxa with abundances 
greater than 1%. In the data set including Emiliania huxleyi, 
only four taxa meet that criterion (Calcidiscus leptoporus, E. 
huxleyi, Gephyrocapsa muellerae and Gephyrocapsa spp. 
small) (Tab. 4.1); whereas in the data set excluding E£. huxleyi, 
seven taxa had abundances greater than 1% (C. /eptoporus, G. 
muellerae, G. oceanica, Gephyrocapsa spp. small, Helicosphaera 
carteri, H. princeiand Umbellosphaera tenuis) (Tab. 4.2). 

The dendrogram based on the data set including EF. 
huxleyi shows two well-defined groups (A and B) distant to 
a compositional Euclidian distance of 12.8 (Fig. 4.1). Group 
A, composed of nine samples (sites 4—6 and 8-13), lies in 
the southeast area and is dominated by G. muellerae (> 56%). 
Group B, composed of 24 samples, is distributed in two 
geographical areas: in the southwest and in the north, and 
is dominated by E. huxleyi (> 93%) (Fig. 5.1). 

The dendrogram based on the data set excluding EF. 
huxleyi, shows two groups (A’ and B’) separated with a 
compositional Euclidean distance of 15.8 (Fig. 4.2). Group 
A’, composed of 14 samples, spans through the southern 


sites and one northern site (sites 1-13 and 14, respectively). 
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TABLE 4 — Relative abundances of the taxa used in the Cluster Analysis. 1) Data set including Emiliania huxleyi. 2) Data set excluding Emiliania 
huxleyi. 
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Figure 4. Dendrograms obtained from the Cluster Analysis (Ward method adapted to compositional data). 1, Data set including £ huxleyi. 2, Data 


set excluding E. huxleyi. 


Group B', composed of 19 samples, corresponds to the 
northern sites. The most abundant species in these groups 
are C. leptoporus and G. muellerae, showing no preference for 
one or the other group. Both are dominant in some sites of 
the group A’ and in some sites of the group B’ (Fig. 5.2). 


DISCUSSION 
Coccolithophore preservation 

Previous to compare the two counting methods, is 
essential to identify if coccolith assemblages are affected 
by dissolution. Therefore, some qualitative and quantitative 
indices have been made. In all the samples, observations in 
light microscope reveal well-preserved coccoliths. The data 
concerning coccolith preservation from different samples 
are given in Table 2. 
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High CEX' index values, all over 0.6, indicate that 
dissolution is not significant. Nevertheless, it is important 
to be careful with the results of this dissolution index, which 
is based on differential dissolution, because it can be 
influenced by the different ecological affinities of the 
species. In the ACM, coccolithophore assemblages tend to 
be E. huxleyi-monospecific northward and this could be 
linked to environmental conditions that influenced the 
ecological preferences of coccoliths. In this area, the CEX’ 
index served as a mere approximation, suggesting an 
integrated approach combining different dissolution proxies 
is necessary. However, Bickert & Wefer (1996) suggested 
that present-day lysocline of calcite in the South Atlantic 
Ocean is at ~ 4,000 meters. All samples included in this 
study were above that depth (Tab. 2). 
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Figure 5. Maps of the Argentinean Continental Margin showing pie-charts depicting coccolith taxa and spatial distribution of assemblages. 1, 
Data set including E. huxleyi. 2, Data set excluding E. huxleyi. A and A’ in black sites numbers; B and B’ in white sites numbers. Pie-charts A, B 
and A’, B’ are the center of assemblages, the numbers indicate the c/r variances. Plotted with GMRT data (Ryan et a/, 2009). 


Coccolithophore assemblages and Cluster Analyses 
Emiliania huxleyi is the only species present in all 
assemblages and, in most of them, is also the dominant 
one. Its highest abundance occurs in the north and in the 
southwest area. Emiliania huxleyi is adapted to a wide 
range of environmental conditions and has a cosmopolitan 
distribution (Winter & Siesser, 1994). This species occurs 
under a broad salinity and temperature range, although its 
highest abundances are at water temperatures of < 23 °C 
(Winter et al, 1979; Haidar & Thierstein, 2001). Brand 
(1994) reported high abundances of E. hux/ey/ in nutrient- 


rich environments such as subtropical central gyres edges, 


upwelling regions and along continental margins. Boeckel 
et al. (2006) observed that £. huxleyi exhibits highest 
abundances in sediments from Brazil, Uruguay and the 
Argentinean Continental Margin (> 80%). Indeed, present day 
records of £. huxleyi blooms have been studied along the 
ACM by Poulton et al. (2011, 2013) and Balch et al. (2014). 
These intense and regular high-reflectance blooms have 
been observed during the austral summer along the ACM 
shelf break and slope, and are enhanced to permanent high- 
chlorophyll zone (Romero et a/., 2006; RUhlemann & Butzin, 
2006; Matano et a/., 2010; Poulton et a/., 2011, 2013). 


If E. huxleyi is excluded from the data set, C. /eptoporus 


73 


Yo § 
APA | Publicacion Electronica - 2023 - Volumen 23(1): 61-77 


(mean 46.45%) and G. muellerae (mean 34.82%) become the 
most abundant species, each dominating different oceano- 
graphic areas: C. /eptoporus is abundant at the southwestern 
sites and at the northern sites, whereas G. muellerae 
dominates the assemblages at the southeastern sites. 
Calcidiscus leptoporus is a cosmopolitan species but is a 
typical component of the cold-water floras (Mcintyre & 
Bé, 1967; Mostajo, 1986). Opinions about nutrient re- 
quirements are diverse, and possibly it could be related to 
the presence of different morphotypes (Knappertsbusch 
et al, 1997), which have not been considered, as it is out of 
the scope of this study. In the southeastern Atlantic Ocean, 
C. leptoporus was associated with intermediate and low 
nutrient conditions in surface waters (Boeckel & Baumann, 
2004). 

Regarding G. muellerae, Holligan et a/ (2010) and 
Charalampopoulou et al. (2016), reported that this species 
is characteristic in the south of the Malvinas Islands and 
shows a preference for very cold and relatively high 
productivity water masses. Likewise, Gephyrocapsa spp. 
small delineates the same trend of G. muellerae, being 
abundant under elevated nutrient conditions (Okada & 
McIntyre, 1977; Saavedra-Pellitero et a/,, 2010; Martinez- 
Sanchez et al, 2019). 

Several environmental parameters are regulating the 
complex ecosystem in which coccolithophores live. Some 
variables of the photic zone, such as temperature, salinity, 
water column stratification and nutrient concentrations, can 
be related to the underlying surface sediment assemblages. 
In other words, coccolith sediment assemblages represent 
the average sea surface environmental conditions of the 
last few hundred years in the area. As more stenotopic 
species are represented in those assemblages, a more 
accurate interpretation of the environmental conditions 
will be possible. Emiliania huxleyi proved to have high 
tolerance to a wide range of environmental conditions and 
is extremely abundant in the studied area. The exclusion of 
this species from statistical analysis might allow to better 
identify geographical distribution patterns of less repre- 
sented, but stenotopic species; or as the environmental 
parameters build up a very complex and synergetic 
ecosystem, the exclusion of £. huxleyi might also represent 


a loss of information. 
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The studied area is restricted to the Subantarctic water 
mass (sensu Mostajo, 1986), although sea surface tem- 
perature difference between the northern and southern 
sites is about 12 °C. In the northern portion of the Drake 
Passage it is about 4 °C and increases northward up to 16 
°C, where the currents separate from the shelf break. How- 
ever, the distribution of the coccolithophore assemblages 
seems to be controlled by other contrasting oceanographic 
features such as nutrient availability or salinity, and not by 
the temperature gradient. Satellite chlorophyll-a observa- 
tions indicate that the shelf break and slope have a high 
nutrient concentration (Romero et a/, 2006; Franco et al, 
2008). Numerical and analytical models suggest that these 
high concentrations are associated with interactions be- 
tween the Malvinas Current and the bottom topography, 
which creates upwelling systems (Piola & Matano, 2017). 

The dendrogram obtained from the Cluster Analysis 
including E. huxleyishows two groups (A and B) and revealed 
a good constraint between the occurrence and distribution 
of species and oceanographic features. Group A, which is 
dominated by G. muellerae, is confined to the southeast area. 
Piola & Matano (2017) suggest that most of the water of 
ACC loops northward to form the MC and a portion deflects 
northward west off the Burdwood Bank. In this area, where 
the MC curve to the east, samples are dominated by G. 
muellerae (group A). As mentioned above, G. muellerae shows 
a good correlation with cool and nutrient-rich environments 
(Boeckel et al, 2006). Additionally, Matano & Palma (2008) 
and Miller et a/ (2011) postulated and verified, using 
simplified models, that the upwelling might be associated 
with frictionally driven intrusions of the MC onto the shelf. 
Matano et al. (2010) observed that this area, where the 
samples are dominated by G. muellerae, is one of the regions 
of largest cross-shelf exchanges. Accordingly, bathymetry 
data revealed that these sites are located near the 200 m 
isobath, where the MC loops. In consequence, elevated 
percentages of G. muellerae could be attributed to nutrient- 
rich waters in this active upwelling region. Group B, which 
is dominated by E. Aux/eyi, includes the southwestern and 
northern sites of the ACM. Emiliania huxleyi, an opportunistic 
species, is also abundant in upwelling areas (Brand, 1994). 
Accordingly, it seems that these two groups could be 
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along nutrient-rich surface waters, where G. muellerae 
exhibit a better response. 

The dendrogram obtained from the Cluster Analysis 
excluding £. huxleyi shows two groups (A’ and B’) which 
cannot be explained by the dominance of the species that 
compose the assemblages at these sites, neither by the 
oceanographic regimen. 

In the ACM, the distribution of the coccolithophore 
assemblages seems to be related to surface nutrient 
availability and this parameter is probably reflected by the 
ecological competition between the most abundant species, 
E. huxleyi and G. muellerae. Further analysis is required to 
better understand the relation between the oceanographic 
regime in the ACM and coccolith assemblages, but the 
exclusion of £. huxleyi from these analyses will incur in an 
important loss of information. 


CONCLUSIONS 

We tested which of the most used data sets showed 
clusters that could be explained by known oceanographic 
conditions in the ACM. The data set including £. huxleyi 
formed a grouping that is geographically well defined and 
could be explained by enhanced upwelling in the ACM, 
whereas the data set excluding £. hux/eyishowed a grouping 
not entirely geographically defined, that could be barely 
explained by the sea surface temperature gradient. In 
addition, the remarkably low number of other species in the 
ACM, in relation to £. hux/eyiand G. muellerae, could limit the 
use of the commonly applied counting methods, excluding 
E. huxleyi. According to the known oceanographic regime 
and models for the ACM, the most appropriate method for 
paleoecological and paleoceanographic reconstructions 
based on recent coccoliths is to consider the assemblage 
including E. huxleyi. This comparison of counting methods 
is a helpful prerequisite for a better understanding of the 
present-day environment influences on the coccoliths 
spatial distribution in the MCA, that will allow us to 
this information in 


extrapolate paleoceanographic 


reconstructions. 
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